Anisotropy-induced ordering in the quantum J\ — J 2 antiferromagnet 



o 
o 



Tommaso Roscilde, 1 Adrian Feiguin, 2 Alexander L. Chernyshev, 2 Shiu Liu, 2 and Stephan Haas 1 

1 Department of Physics and Astronomy, University of Southern California, Los Angeles, CA 90089-0484 
2 Department of Physics and Astronomy, University of California, Irvine, CA 92697 

(Dated: February 2, 2008) 



We study the effect of spin anisotropics on a frustrated quantum antiferromagnet using the Ji- 
J xxz model on the square lattice. The T = and finite- T phase diagrams of this model are 
obtained utilizing spin-wave theory, exact diagonalization, and quantum Monte Carlo. We find that 
anisotropic frustration tends to stabilize XY- and Ising-like ordered phases, while the disordered 
spin-liquid phase is restricted to a small region of the phase diagram. The ordered phases are 
separated by first-order transitions and exhibit a non-trivial reentrance phenomenon. 
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Frustrated quantum magnets attract significant inter- 
est because of the spin-liquid and novel ordered phases 
they may exhibit. Q, 0, The macroscopic degeneracy 
of the ground state in such magnets makes them very 
sensitive to additional interactions that may lead to var- 
ious unconventional ordered states. Perhaps the most 
extensively investigated, yet highly controversial model 
is the 5=1/2 J\-Ji Heisenberg model on the square 
lattice with competing nearest-neighbor, J\, and next- 
nearest- neighbor, J2, interactions. |4j The recent synthe- 
sis of compounds that can be closely described by the 
two-dimensional (2D) J\-Ji Hamiltonianjsj has also fu- 
eled interest in the properties of this basic model of mag- 
netic frustration. In that respect, the presence of spin 
anisotropics in real systems raises the question of how ro- 
bust the behavior of the isotropic J1-J2 model is against 
such perturbations. 

In the present paper we study the effect of spin 
anisotropies on the properties of frustrated quantum an- 
tiferromagnets using a generalization of the J\-Ji model, 
in which the frustrating next-nearest-neighbor interac- 
tion is anisotropic. Such a J\-J xxz model is given by: 



E 



jz nz nz 
J 2 D i D k ' 



■J 2 S i 



(1) 



where Si is the spin operator, the sites k are on the 
square lattice, and the summation runs over all nearest- 
neighbor (n.n.) or next-nearest-neighbor (n.n.n.) sites. 
We hereafter use the dimensionless ratios a z = J%/Ji 
and a± = / Ji ■ 

For the extensively studied isotropic case a± — a z , 
analytical and numerical works have suggested the ex- 
istence of a T — 0, non-magnetic gapped phase for 
0-4 ;$ -h/Ji < 0.6, separating the Neel state from the 
collinear state. |^ The absence of long-range order at any 
finite temperature in 2D systems for the isotropic Ji- 
J2 case is dictated by the Mermin- Wagner theorem. In 
the collinear phase, however, the 2D J1-J2 model has 
been predicted to exhibit an Ising-like Chandra-Colcman- 
Larkin (CCL) transition Q at T > with spontaneous 
breaking of the discrete lattice rotational symmetry. This 
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FIG. 1: Zero-temperature phase diagram of the 5 = 1/2 
Ji-J xxz model within linear spin-wave theory. Solid and 
dashed lines denote the second- and first-order transitions, 
respectively. Inset: non-magnetic region magnified. 



unusual scenario has recently received support from nu- 
merical and analytical studies. 

In this Letter we show that even a small anisotropy 
in the frustrating (n.n.n.) coupling leads to strong de- 
viations from the behavior of the isotropic J1-J2 model. 
The main effect of anisotropy is to energetically favor 
the sector of spin space in which the frustration is min- 
imal. This tendency of spins to avoid frustration leads 
to somewhat counter-intuitive results, as, e.g., Ising-like 
anisotropy in J2 favoring the XY phase and vice versa. In 
the case of the unfrustrated 2D Heisenberg antiferromag- 
net, an arbitrarily small anisotropy is found to stabilize 
finite-temperature ordered phases. |2j In this work we 
show that the anisotropically frustrated J\-J x model 
is generally more ordered than its isotropic counterpart. 

The T = phase diagram of classical (S = 00) 
Ji-J xxz model is shown in Fig. ^ We find four 
phases characterized by different order parameters, = 
(%2 r e, lqr S"), hereafter indicated as: (i) XY-Neel phase 
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rri?^ j, (ii) Z-collinear phase with mf^ 
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(Hi) XY-collinear phase with rnf% m or m^^, and (iv) 
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Z-Neel phase with m^. ^ . The diagonal line a± = a z 
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corresponds to the isotropic J\-J% model. All four phases 
are separated by first-order transition lines. Taking into 
account quantum corrections within linear spin-wave the- 
ory results in the opening of a spin-liquid region where 
all the order parameters vanish. Along the main diagonal 
this reproduces the well-known results for the isotropic 
case.y Our analysis suggests that in the vicinity of the 
a±_ = a z line the T = transitions from the ordered 
phases to the spin-liquid phase is of second order. How- 
ever, we also find that the collinear phases are separated 
from the spin- liquid phase mostly by the first- order tran- 
sition lines, as shown in Fig. ^ This means that the 
XY- and Z-collinear orders are particularly robust, so 
that even quantum-renormalized order parameters van- 
ish with a finite jump when crossing such boundaries. 

Given the anisotropic nature of the ordered phases 
found in the T = phase diagram, one can anticipate 
qualitatively the finite temperature behavior. In the Ising 
regions magnetic order is expected to survive at T > 
up to a second-order transition point. Topological order 
is instead expected at T > in the XY regions up to a 
Berezinskii-Kosterlitz-Thouless (BKT) transition. 

The study of the quantum J^-J** 7, model at finite 
temperature by quantum Monte Carlo (QMC) techniques 
is generally precluded by the sign problem [2j due to the 
transverse frustrating coupling a±. However, by remov- 
ing that term completely we obtain the sign-problem- free 
J1-J2 model, where the frustration is Ising-like, corre- 
sponding to the a± — line of the phase diagram in Fig. 
n A comprehensive study of this limiting case is, there- 
fore, very important since its thermodynamics can shed 
light on the main features of a large region of the phase 
diagram (a± < a z ). Moreover, in this limit of Jf-only 
frustration, the Hamiltonian of Eq. Q is also relevant 
for the broad class of strongly correlated bosonic systems. 
In fact, this J\-J% model can be mapped exactly onto the 
frustrated hard-core Bose-Hubbard model. 

We have thus investigated the thermodynamics of 
the S — 1/2 J1-J2 model by means of the Stochas- 
tic Series Exp ansion QMC based on the directed-loop 
algorithm. JlOj A large interval of frustration values, < 
a z < 1.7, has been scanned using different lattices L x L 
up to L = 96 to perform an extensive finite-size scal- 
ing analysis. The resulting phase diagram of the J\-J% 
model is shown in Fig. [21 As main features we observe 
the occurrence of the topologically ordered phase (quasi- 
long-range XY"-Neel) for small a z and the zT-collinear or- 
der for large a z . These two phases are separated by a 
first-order transition line. Phases with finite-temperature 
topological order and zT-collinear order in the spin lan- 
guage correspond, in the bosonic language, to quasi-long- 
range superfluidity and striped solid order, respectively. 

The study of the first-order transition line has been 
performed by making use of a quantum parallel temper- 
ing technique |ll| to overcome critical slowing down at 
first-order transition points. 01 I n the collinear phase 
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FIG. 2: Phase diagram of the S = 1/2 Ji-Jf model. Mag- 
netic (and corresponding bosonic) phases are indicated. The 
horizontal dashed line marks the BKT critical temperature of 
the .5 = 1/2 XY model |l5fl . the dotted line is a logarithmic fit 
for small a z (see text). Inset: zoom on the tricritical point. 

we have also used thermal parallel tempering to over- 
come the well-known loss of efficiency of standard cluster 
algorithms in Ising-like frustrated antiferromagnets. [LH 

A remarkable feature of the phase diagram is that even 
an infinitesimal frustration, a 2 « 1, induces a finite- 
temperature BKT transition. The critical temperature 
Tbkt has been extracted through scaling analysis of the 
helicity modulus T and of the transverse structure factor 
gxxiyy)^ y 0, E| As expected for weakly anisotropic 
systems, jjj Tbkt depends logarithmically on the frus- 
tration as Tbkt ~ 47T/9 S /(C + |lna z |), where p s is the 
spin stiffness of the unfrustrated Heisenberg antiferro- 
magnet and C is a constant. A fit of our data to the 
above law for a z < 0.1 gives p s = 0.175(2) in excel- 
lent agreement with the best available estimates, [l6| 
and C = 5.4(1). For larger frustration, the BKT transi- 
tion temperature reaches a maximum which falls very 
close to the critical temperature of the S — 1/2 XY 
model, T BKT /J = 0.3427(2), indicating that the J x - 
J2 model with frustration a z ss 0.65 is an almost ideal 
realization of the quantum XY model. 

At zero temperature, the XY-Neel and the Z-collincar 
phases are separated by a first-order transition at the 
critical value a* = 1.252(5). We notice that quantum 
fluctuations cause a 25% shift of this critical value with 

(c) 

respect to the classical value a\ ' = 1, thus promoting the 
Neel phase against the collinear phase. When increasing 
the temperature the first order transition line shows a 
reentrant behavior, as seen in the inset of Fig. [5] For 
1.21 < a z < 1.25, when temperature is increased, the 
topological order becomes unstable to the onset of the 
collinear magnetic order, which carries a higher entropy 
content. One can argue that the bending of the first-order 
transition line reflects the reduced role of quantum flue- 
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FIG. 3: Left panel: the double-peak structure in the his- 
tograms of the collinear magnetization m co ii for system size 
L — 32 and temperature T/Ji = 0.36. Right panel: Binder's 
fourth cumulant of the collinear magnetization and the CCL 
order parameter for a z = 1.3. 

(c) 

tuations, responsible for shifting the critical value a z , 
as the temperature is raised above zero. We notice that 
similar phase reentrance phenomena have recently been 
observed in related magnetic and bosonic systems. [l7| 

The first-order transition line meets the BKT line at 
the point (a z ,T/Ji) = (1.210(5), 0.28(1)), above which 
it maintains its first-order nature up to a temperature 
T I J\ k, 0.36. This is shown in Fig. 03 where the phase 
coexistence manifests itself in the double-peak struc- 
ture of the distribution for the collinear magnetization 
"i-coii = mfff o) + m fo tt) aroun d the transition point. We 
have checked that the double-peak feature persists when 
increasing the lattice size up to L = 48. For T > 0.36Ji, 
determining the order of the transition to collinear order 
becomes more complicated. In fact the distribution of the 
order parameter m co \\ loses the strong two-peak feature, 
suggesting the transition to be weakly first- order. For 
a z ^> 1 the model reduces to two disconnected Ising mod- 
els, living on the two sublattices of the original square 
lattice. In this limit, according to Onsager's solution, 
we should expect a 2D Ising transition at a critical tem- 
perature T co ii = 2JfS'7ln(l + y/2) = 0.567 Jf. The 
crossover from a weakly first-order to a second-order 2D 
Ising transition might happen with a continuous change 
of the critical exponents for increasing a z , as argued to 
occur in the frustrated 2D Ising antiferromagnet. 

We have determined the transition points T co \\ for T > 
0.36 by the crossing of Binder's fourth cumulant = 
1 — (m* oll ) /3(ml oll ) 2 . In Fig. ||we observe that, for 
a z — 1.3 as well as for the other cases studied, the values 
of U^' (L) for different lattice sizes do not cross at the 
2D Ising critical value C 4|C = 0.6107 0, which suggests 
the universality class of the transition to be distinct from 
2D Ising. 

An interesting question is whether the CCL transition 
of the isotropic J1-J2 antiferromagnet, with sponta- 
neous breaking of the discrete lattice rotational symme- 
try, survives in the collinear phase of the anisotropic 
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FIG. 4: Zero-temperature phase diagram of the Ji-J xxz 
model from exact diagonalization on clusters of size N. 

J1-J2 model. The order parameter for the transition 
reads: i @ = - S? +ltj+1 )(Sf +1<j - 5f J+1 ))/4. 

We have studied the scaling of Binder's fourth cumu- 
lant related to the order parameter er, U± = 1 — 
("coll) /3( <T coii) 2 - A s shown by the crossing of the Binder's 
cumulant in Fig. |3 the CCL transition seems to occur 
at the same (or slightly higher) temperature with respect 
to the magnetic transition. 

Let us now return to the general Ji-J** 2 model 
to investigate how the ground state properties evolve 
in presence of the transverse component of frustration, 
Q!_l < ol z . Making use of exact diagonalization on lat- 
tices of iV =16, 20, and 32 spins we have determined the 
dependence of the order parameters, mj? j and m co u, on 
a z at fixed a±. Using a dense mesh of points, we obtain 
a qualitative finite-size estimate of the phase boundaries 
from the vanishing of the second derivatives c^m"^/ 'da 2 , 
marking an inflection point in the order parameters. In 
this way we have obtained the T = phase diagram 
shown in Fig. 0] For sufficiently low transverse frustra- 
tion (aj_ < 0.4) both order parameters display an inflec- 
tion at the same point, suggesting that the transverse 
frustration does not introduce any intermediate phase. 
The sharpness of the peaks in the derivatives of the or- 
der parameters suggests that the transition maintains the 
first-order nature found in the Ji-JJ model. This first- 
order transition line runs roughly parallel to the classical 
line, being shifted to higher a z by about 20% due to 
quantum fluctuations. As the transverse frustration is 
increased, the inflection points of the order parameters 
separate from each other, leaving space for an interme- 
diate region with no magnetic order, continuously con- 
nected with the spin- liquid phase in the isotropic limit. 
Nonetheless, we observe the intermediate non-magnetic 
phase to shrink dramatically as the cluster size is in- 
creased. The small size of the clusters does not allow 
us to conclude on the real extension of the intermedi- 
ate non-magnetic phase and on the nature of the phase 
transition line(s) for a± > 0.4. Further investigations of 
this issue are currently in progress. A suggestive scenario 
is that the transition line between the collinear and the 
non-magnetic phase remains of first order all the way to 
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FIG. 5: Global phase diagram of the J x - J* xz model. The 
leftmost filled area represents the region with broken rota- 
tional symmetry in lattice space but unbroken symmetry in 
spin space; data for the CCL transition are from Ref. 0. 

the isotropic li mit, in agreement with the predictions of 
previous works 20]. 

We finally propose a global phase diagram for the 
Ji-J* xz model, shown in Fig. [SJ by merging the 
the T — data with the finite-temperature data. 
Here the finite-temperature Z-collinear and topological 
phases are separated by a first-order transition surface, 
and upper-bounded by a BKT transition surface and 
a weakly- first-order/second-order surface, respectively. 
For a z = a± the above critical surfaces vanish. More- 
over, in the high- frustration region of a± « a z ~ 0.5 
we adopt the picture of an extended spin-liquid region 
where all transition temperatures are also vanishing. 
For a± — > a z outside the spin-liquid region, the crit- 
ical temperatures T c vanish with a very steep slope, 
following a logarithmic dependence on the anisotropy 
A = \a z - a±\ < 1, T c ~ llnfA)]- 1 . This depen- 
dence is obtained through the mean- field condition |7ll21| 
ksT c ~ JiA[£( Jl ~ j2 )] 2 , assuming an exponential diver- 
gence of the correlation length of the isotropic Jy-J^ 
model, £( J i-^) ps T _1 exp(27ro^ 1 " J2) /T 1 ) (as found, e.g., 
by modified spin- wave theory |22(). At the spin-liquid re- 
gion boundaries, instead, the spin stiffness p s vanishes; to 
account for that, along the boundary between the spin- 
liquid and the Neel region we introduce a A-dependent 
stiffness in the mean- field equation, p s ~ (A — A c )^, 
where A = A c marks the boundaries of the spin-liquid 
region. In drawing the phase diagram we have taken 
(3 = 1, as suggested by the data in Ref. HjJ In this 
way we get a much smoother vanishing of T c , T c ~ 
(A — A c )/| ln(A — A c )]|. For the isotropic case a z = a± 
and in the collinear phase, the only finite-T transition 
line is the CCL one, recently estimated for S — 1/2 in 
Ref. 0- Having observed the CCL transition to occur 
very close to the magnetic transition for a± = 0, we ex- 



pect the CCL transition surface to quickly merge with 
the steep magnetic transition surface as a± < a z . Thus, 
we suggest that the CCL transition exists as a distinct 
thermodynamic feature only in a very close vicinity of 
the isotropic limit. 

In summary, we have shown that the Heisenberg anti- 
ferromagnet on the square lattice with anisotropic next- 
nearest-neighbor frustration shows a very rich phase di- 
agram characterized by several finite-temperature tran- 
sitions, which we studied using QMC, exact diagonaliza- 
tion, and spin-wave theory. The competition between or- 
dered phases leads to non-trivial reentrance phenomena, 
which is also common to related bosonic systems. Since 
anisotropic spin-spin interactions are a general feature of 
real magnets, the ordering effect due to anisotropic frus- 
tration is a realistic mechanism to explain magnetic tran- 
sitions in recently synthetized frustrated quantum anti- 
ferromagnets. 
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